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Abstract 

We propose an (essentially combinatorial) approach to the correlation 
functions of the domain wall six vertex model. 

We reproduce the boundary 1-point function determinant expression 
of Bogoliubov, Pronko and Zvonarev, then use that as a building block to 
obtain analogous expressions for boundary 2-point functions. 

The latter can be used, at least in principle, to express more general 
boundary (and bulk) correlation functions as sums over (products of) 
determinants. 



1 Introduction 

Computing off critical correlation functions 1 is probably the most challenging 
open problem currently under investigation in exactly solvable lattice models 

m 

The six vertex model, with domain wall boundary conditions (dwbc's) is an 
ideal testing ground of possible approaches to such computations, particularly 
if one is interested in computations on a finite lattice 2 . 

The model was first introduced by Korepin [2], who also formulated recursion 
relations that uniquely determine the partition function. Korepin's recursion 
relations were solved by Izergin [3], who obtained a determinant representation 
for the partition function. 

'Supported by the Australian Research Council (ARC). 

1 The literature on off critical correlation functions is extensive. We refer the reader to [1] 
for references to the literature up to the early 90's and to a search of |http : //arXiv . org| for 
more recent literature. 

2 The previous footnote applies verbatim to the literature on the domain wall six vertex 
model. 
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Using the algebraic Bethe ansatz py> Bogoliubov, Pronko and Zvonarev in- 
troduced a definition of boundary 1-point functions in the model |o], and com- 
puted them in determinant form 3 . 

Inspired by Bogoliubov et al, we propose a definition of boundary n-point 
functions, that in the specific case that they considered, (almost) coincides with 
theirs (up to obvious factors due to the differences in our definitions of what a 
boundary 1-point function is). 

We use combinatorics, based on the Yang Baxter equation, to compute deter- 
minant representations for boundary 1-point, and 2-point functions, and outline 
the (almost mechanical) extension to (certain) boundary n-point and bulk cor- 
relation functions. In the case of boundary 1-point functions, we reproduce the 
result of Bogoliubov et al [5] using elementary manipulations. 

Our derivations are (basically) combinatorial and rely on repeated applica- 
tion of graphical operations, hence the proliferation of figures in the paper 4 . 

There has been interest in this very model from algebraic combinatorialists 
over the past decade, particularly since the work of Kuperberg . Our exposi- 
tion is elementary in the hope that it will serve as an introduction to this part 
of exactly solved statistical mechanical models for non-physicists. 

2 The model 

To be reasonably self-contained, we start by recalling basic definitions related 
to lattice models in general, and to the dwbc six vertex model in particular. We 
refer the reader to ^ 0] for further details. 

Lattice configurations Consider a square lattice, with N c vertical lines 
(columns) labeled from left to right as {N c , N c — 1, • • • , 2, 1}, and N r horizontal 
lines (rows) labeled from top to bottom as {1, 2, • • • , N r — 1, N r }. Initially, we 
take N c — N r = N. Later, we will relax this condition and consider lattices 
with deviations from dwbe's and N c ^ N r . 

1 

2 

3 



N 

N 3 2 1 
Figure 1 Our playing field 

'^Bogoliubov et al obtained determinant expressions for the 1-point functions, and also for 
the boundary spontaneous polarization, which is a more difficult problem. Here, we discuss 
only the former. 

4 Our results can be produced using the algebraic Bethe ansatz. However, our method is 
elementary, and hopefully there is virtue in having more than one approach to this problem. 
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Boundary lines or simply 'boundaries', are the right most, left most, top 
and bottom lattice lines. 



Orientations We assign each lattice line an orientation. In our convention, 
horizontal lines are oriented from left to right. Vertical lines are oriented from 
bottom to top. 



Figure 2 The white arrows denote the orientations of the lattice 
lines. 



Rapidities We assign each oriented lattice line a complex variable called a 
rapidity. We assign the i-th horizontal line a rapidity x i} and the j-th vertical 
line a rapidity yj . We restrict our attention to the fully inhomogeneous situation 
in which all rapidities are different. 



Xi 

x 2 



xn- 



flHI 



Vn 



2/2 2/i 



Figure 3 The rapidity flows along the lattice lines. The direction 
of the flow is the orientation of the line. The strength of the flow is 
the rapidity. In our convention, the flow in each line orthogonally 
crosses lines with decreasing label. 

Bonds are horizontal or vertical line segments that lies between two intersec- 
tion points, or at the very end of a line. 
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Arrows To each bond we assign an arrow that can point in either direction 
along the bond. 

Vertices We call the intersection point of the z-th horizontal line and the j-th 
vertical line, together with the 4 bonds attached to it, and the arrows on these 
bonds, a vertex My. 

This extends to the general case when the two intersecting lines are not 
'everywhere' straight. The line orientations are sufficient to determine which of 
the two intersecting lines is 'locally' vertical and the other is 'locally' horizontal 5 . 

Weights We assign each vertex, i>y , a (Boltzmann) weight, Wy, that depends 
on the following data 

1. The orientations of the lines that intersect at that vertex. If necessary, we 
must rotate a vertex so that one line is oriented from left to right, while 
the other is oriented from bottom to top. 

2. The orientations of the arrows on the four bonds attached to the intersec- 
tion point. 

3. The rapidities, {xi,yj}, flowing through the intersection point. 

4. A parameter A, called the crossing parameter that is the same for all 
vertices Vij , and characterises the physical characteristics of the model. 

The six vertex model Given that each vertex has 4 bonds, and each bond 
carries an arrow with 2 possible orientations, there are 16 possible vertex con- 
figurations. The six vertex model is the special case where only six types of 
vertices generally have non-zero weights. The weights of the remaining 10 types 
are set to zero, weights. 




t^y §y A y 

a(x,y) b(x,y) c(x,y) 



Figure 4 The vertices of the six vertex model, with their rapidity flows. 
The weights of every two vertices in the same column are equal, and is 
shown below it. 

5 When in doubt, deform the intersecting lines locally such that their orientations agree 
with those given in figure HI 
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Conservation of arrow flow The vertices with non-zero weights are pre- 
cisely those that conserve arrow flow: In each vertex, two arrows point inwards, 
and two points outwards. 

Bracket notation We define [x] = sinh(Ax), where A is the crossing param- 
eter. 

Vertex weights In the bracket notation, the vertex weights are 



a{xi,Vj) = \-Xi + yj + 1], b(xi,yj) 



+ Vj], c(xi,Vj) = [l] (1) 



where Xi (respectively yi) is the horizontal (respectively vertical) rapidity vari- 
able flowing through the vertex. 

The weights of type a and type b vertices depend on the differences of the 
rapidities flowing through the vertex, and can vanish for specific values of the 
rapidities. The weight of the c vertex is independent of the rapidities. 

Yang Baxter equations The origin of solvability of the six vertex model 
(irrespective of the boundary conditions) is that the (Boltzmann) weights of the 
model satisfy the Yang Baxter equations pQ. Defining the (R) matrix 



a(x,y) 

b(x,y) c(x,y) 

c(x,y) b(x,y) 





R{x,y) = 

the set of all Yang Baxter equations can be written in matrix form as 







a(x,y) 



(2) 



R(x, y)R(x, z)R(y, z) = R(y, z)R(x, z)R(x, y) (3) 

where matrix multiplication is implied. The matrix equation |21 is equivalent to 
a set of equations between the weights. An example of a Yang Baxter equation, 
that will be used below, is shown graphically in figure El 




+ 









A 


X 


y 


z 


X 





M ♦ M ♦ M ♦ 

Figure 5 A Yang Baxter equation. The external arrows, rapidities, and 
line orientations are the same on corresponding lines. Line orientations 
determine the vertex type unambiguously. 
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The Yang Baxter equation corresponding to figure [3] is 



b(y, z)a(x, z)c(x, y) + c(y, z)c(x, z)b(x, y) = c(x, y)b{x, z)a(y, z) (4) 

Domain wall boundary conditions (dwbc's) For a vertex model on a 
finite square lattice, with N c = N r = N, one can require that all arrows on the 
left and right boundaries point inwards, and all arrows on the upper and lower 
boundaries point outwards. These are called 'domain wall boundary conditions' 
(dwbc's) 6 . An example is shown in figure 




Partition function Given a statistical mechanical model, the partition func- 
tion is defined as the sum over all weighted configurations. The weight of each 
configuration is the product of the weights, Wij, of its vertices My 

z N ({x},{y}} = E II w a ( 5 ) 

all allowed vertices 
configurations 

Remark 1 The expression in equation is computationally worthless, as it 
cannot be evaluated in polynomial time in N . 



2.1 Izergin's expression for the partition function 

In Izergin solved Korepin's recursion relations 0, and obtained an explicit 
expression for the partition function of the dwbc six vertex model in terms of a 
determinant 

Zn {{x},{y}j = n b(x . )b( . detM 



6 One can equally well make the opposite choice of boundary arrow orientations. For 
consistency, we need to be make one choice and stay with that. This work is full of such 
choices. 
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nSf=i[-«* + vj + /«x 

= Tpj-^ = T7 7 detM (6) 

where 

M = c{xi,yj) = [1] 

lJ a(xi,yj)b(xi,yj) [-Xi+yj + l][-x i + y j \ 

Remark 2 Izergin's expression, equation^ for the partition function in terms 
of a determinant is computationally meaningful, as it can be evaluated in poly- 
nomial time in N . 



3 Correlation functions 

Apart from the partition function which is the weighted sum over all configura- 
tions allowed by the dwbc's, one wishes to consider n-point correlation functions. 
These are defined as weighted sums over all configurations such that n arrows 
(that are summed over in the partition function) have certain frozen orienta- 
tions. 

Since the arrows on the outside of all vertical and horizontal lines are always 
frozen by definition of dwbc's, we need to consider freezing arrows that are on 
the inside of the lattice. 

We start by considering the 1-point function obtained by summing all con- 
figurations such that a certain single arrow is kept frozen. Without loss of 
generality, let us consider freezing a horizontal arrow. 

For example, consider freezing the orientation of the white arrow in figure 
where bonds without arrows indicate bonds whose arrow orientations are 
summed over 



1 - 

2 - 

3 - 

N- 



N 



1 



Figure 7 The configurations correspond- 
ing to a 1-point function. 

Our proposal is that one should start by slicing all such configuration into 
two parts, a right part, and a left part 
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Figure 8 Slicing a domain wall configuration into two parts. The 
white arrows are frozen. The shaded exposed bonds indicate iden- 
tified arrows that are summed over. 

where the dashed lines in figure |H1 stand for bonds whose arrows are not frozen, 
and can take any orientation that is allowed. 



Exposed arrows We call the arrows that were originally internal, and be- 
came external upon slicing the lattice, 'exposed' arrows. When we slice an 
NxN lattice into two parts, there will be N exposed arrows on each part. 
The orientations of two exposed arrows that originate from slicing an originally 
internal bond are identical, and (generally) summed over. 



Boundary correlation functions We wish to define each part, obtained 
by slicing a set of domain wall configurations, as in figure [SI as a boundary 
correlation function, with the understanding that arrows that were frozen by 
the original dwbc's remain as they were, n exposed arrows (on each part) are 
frozen, the rest of the exposed arrows are summed over. In the example in figure 
El n= 1. 



Simple boundary n-point functions There are special cases when apart 
from the n frozen arrows, the rest of the arrows assume domain wall bound- 
ary orientations: all those on the outside of the vertical boundary lines point 
inwards, and all those on the outside of the horizontal boundary lines point 
inwards. We refer to these as simple n-point boundary functions. For example, 
each part in the following set of configurations (with all internal arrows summed 
over) is a simple boundary 2-point function: 
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N 3 2 1 

Figure 9 Two parts, each of which is a simple boundary 2-point 
function. 

Origin of the simple boundary functions From dwbc's, and conservation 
of arrow flow, there is exactly one horizontal arrow, pointing to the right, on the 
left of the right boundary, two such arrows to left of the second column, three 
to the left of the third, and so forth, until we reach the left boundary and find 
the correct dwbc's. 




N 3 2 1 



Figure 10 Scanning the white arrows from right to left, we can see 
how the number of horizontal arrows pointing to the right increases 
by 1, every time we step to the left. 

Simple n-point functions, with n horizontal arrows, all on the same column, 
are obtained by slicing a domain wall configuration vertically, between the n-th 
and vertical line (counting from the right) and the (n+l)-st line. 

Precisely the same arguments apply to horizontal arrows that point to the 
left, and vertical arrows that point cither up or down. 



Composite boundary ?i-point functions are obtained, for example, when 
an iVx N lattice is sliced, N r i g ht exposed arrows are pointing to the right, N — 
N r ight are pointing to left, and we freeze n arrows that point, let's say, to the 
right. If n = N r i g ht, then the boundary function is simple. If n < N r i g ht, we 
need to sum over the positions of the remaining right-pointing arrows. Hence, 
such a boundary function is not simple. 

Bulk correlation functions In principle, bulk n-point function, where all 
arrows point in the same direction, and lie on the same line, are obtained by 
taking products of boundary functions corresponding to each part. 

More general bulk functions can be obtained by taking products of more 
general boundary functions. We have nothing deep to say about bulk functions 
at this stage. 

4 Simple boundary 1-point functions 

The simplest correlation function to compute is a simple boundary 1-point func- 
tion. This will correspond to freezing any one of the white arrows shown in figure 

El 




N 3 2 1 



Figure 11 There is exactly one internal arrow, that touches a 
boundary, and that points opposite to the external arrows on the 
other side of that boundary. 

We choose to freeze the horizontal right-oriented arrow. Using the algebraic 
Bethe ansatz, Bogoliubov et al jH] obtained a determinant representation of this 
correlation function. More precisely, they obtained the partition function of 
the set of all domain wall boundary configurations such that the c vertex on 
the right boundary lies on row r, normalized by the full partition function Zjy. 
This is their definition of a boundary 1-point boundary function. 

To reproduce their result, we start by slicing the lattice vertically, into two 
parts, and expose the frozen arrow 
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Figure 12 Slicing a domain wall configuration into two 1-point 
functions. 

We end up with two simple 1-point functions: One on the right, and one on 
the left. The right 1-point function is can be computed by inspection. The left 
will be computed below. The result of Bogoliubov et al is the product of the 
two parts, normalized by the partition function. 



4.1 The right boundary 1-point function 

Consider the partition function of the lxiV right part, and take the single c 
vertex to be on row r. 

Given the dwbc's and conservation of arrow flows, all vertices that are above 
row r, will be of type b, all those below will be of type a, and the weight of the 
right part is 



' J K x r,Vl) 



N 



i=r-\-l 



(7) 



4.2 The left boundary 1-point function 

Consider the partition function of the (N— l)xN left part. If the inverted 
arrow on the right boundary were all the way at the top row, we would have 
the situation shown in figure IT^l 



11 



N 3 2 



Figure 13 'Peeling' a frozen row. 

In this case, all vertices on the top row are frozen to be of type a. Therefore, 
we can peel that top row, and end up with configurations on an (JV— l)x(JV— 1) 
square lattice that satisfy dwbc's, and whose partition function can be evaluated 
using Izergin's expression. 

However, the inverted arrow is, in general, not at the top row. To bring it 
to the top row, we use the Yang Baxter equation as follows. 

4.3 Rolling once 

Consider the set of all configurations that remain after peeling the right bound- 
ary. They live on an (N-l)xN lattice, with dwbc's, except for an inverted 
arrow at row r, on the right boundary. 

Take the inverted arrow to be (initially) at row i, and take the horizontal 
rapidity that flows through it to be (initially) Xi. We denote the partition 
function of the set of all such configurations by F[ri,Xi], and represent them 
schematically by the graph on the left hand side of figure ITU 

Consider another set of configurations, F[ri—±, Xi—i], that is identical to 
F[ri,Xi], except that the inverted arrow is now on row r,-i, and represent it 
schematically as the graph on the right in figure PHI 




Figure 14 The graph on the left is a schematic presentation of 
the set of configurations with an inverted arrow that we wish to roll 
upwards. The figure on the left is an identical set of configurations, 
apart from the fact that the inverted arrow is one row higher. 
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Multiplying the graph on the left by b(xi-\, Xi), and the graph on the right 
by c(xi-i,Xi), and adding the results 




Figure 15 Multiplying the graph on the right by a type b vertex 
and that on the left by a type c vertex and adding the results. 

Since the arrows in the external loop on the right have both possible orien- 
tations, they are summed over, and we have 




Figure 16 All possible arrow configurations that are allowed in the 
external loop on the right are summed over. 

But now we are in a position to use the Yang Baxter equation to move that 
external vertex, horizontally through the lattice, all the way to the left hand 
side 




Figure 17 The result of using the Yang Baxter eguation to weave 
the external loop all the way to the left. All possible arrow config- 
urations in that loop are summed over. 

Given the orientations of the boundary arrows on the left hand side, the 
vertex that emerges is uniquely a type a vertex 
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Figure 18 The only possible external vertex on the left is an a 
vertex. 

Equating the initial sum with the final result, we obtain 

F[ri,Xi]b(xi-i,Xi) +F[ri-i,Xi-i]c(xi-i,Xi) = a(x i -i,Xi)F[r i ^i,x l ] (8) 

Dividing both sides by b(xi—\, xi), moving the second term on the left 
hand side to the right, and using b(xi,Xi—i) = —b(xi—i,Xi), and c(xi,Xi-i) = 
c(xi-i, Xi), we end up with the identity 

F[ri,Xi] =f(xi-i,Xi)F[r i - 1 ,Xi]+g(xi,Xi- 1 )F[ri-i,Xi- 1 ] (9) 
where we have defined 

_ a{xi,x l+ i) _ c(xj,x i+1 ) 

Ji,i+1 — , / \ ? <?M+1 — i / \ 

b(Xi,x l+ i) b(Xi,x i+ i) 

Remark 3 We write the weight of a type c vertex with formal dependence on its 
rapidities to make it easier to keep track of its position, although it is independent 
of the rapidities. 

Equation El says that we can rewrite the set of all configurations, in which 
the inverted arrow is on row i, in terms of configurations, in which the inverted 
arrow is in row (i— 1). 

We have succeeded in 'rolling' the inverted arrow by one row towards the 
top. The price we have to pay for that is that we end up with two sets of 
configurations instead of one. 

Further, it easy to see that, if the initial inverted arrow is at row r, then 
rolling it upwards as we did above, the number of configurations will keep on 
doubling, leading to 2 r ~ 1 configurations. 

However, as we will see below, the Yang Baxter equation can be repeatedly 
used to reduce this number from 2 r ~ 1 to r, which is typical of how the Yang 
Baxter equation makes models solvable 

To simplify the notation, we define 

F%=F[ ni ,x i2 } (10) 

as the configuration with the arrow on the ii-th horizontal line inverted, and 
the rapidity through that line is x%, so we have 

F* = f i . lii F*_ 1 + g^ 1 F*zi (11) 
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4.4 Rolling twice 

Suppose we roll an inverted arrow up twice. From equation^] we obtain 



Ft = fi,i-\F$-l + 9i-l,iFi-i 

— fi,i-l(fi,i-2Fi_2 + 5i-2,i^T_ 2 2 ) 

+ 9i-i,iUi-i,i-2FlZ 2 + 9i-2,i-iF-Z 2 ) 

= fi.i-lfi,i-2F-_ 2 + fi-l t i-2gi-l,iF-_ 2 
+ (0i-2,i/i,i-l + gi^l,i9i~2,i-l)FlZ2 (12) 

Equation 1121 can be simplified by combining the two terms in the coeffi- 
cient of -Fjrf into one using a Yang Baxter equation as follows. Re-expand the 
coefficients of Ffzi m terms of Boltzmann weights to obtain 

9i-2,ifi,i-l + 9i-l,i9i-2,i-l = (13) 

c(xi-2,Xi) a(xi,Xi-i) c(xi_i,x,) c(xi- 2 ,Xi-i) 



b(xi- 2 ,Xi) b(xi,Xi-i) b(xi-i,Xi) b(xi-2,Xi^i) 



(14) 



Rewriting the right hand side of eauation ll3l in a slightly more convenient form, 
we obtain 

(15) 

c(xi,Xi_2) a(xi,Xi-i) c(xi,Xi-i) c(xi-i,Xi-2) 



b(xi,Xi-2) b(xi,Xi-x) b(xi,Xi-i) b(xi-i,Xi-2) 

Multiply both terms in equation ^] by 



(16) 



b(x i: Xi-i)b(xi,Xi-2)b(xi-i,Xi-2) (17) 

we obtain 

- a(xi,Xi-x)c(xi,Xi-2)b(xi-i,Xi-2) + c(xi,Xi-x)c(xi-i,Xi-2)b(xi,Xi-2) (18) 

Comparing eg uat ion 1 1 81 with the left hand side of the Yang Baxter equation 0] 
we find that they are identical, if we identify x — > i,y — *■ £ — 2,z— * i — 1. Using 
these identifications, we can write equation 01 as 

c(xi-l, Xi-2)c(Xi, Xi-i)b(Xi, Xj_ 2 ) + b{Xi-2,X i -{)a(xuX i - 1 )c{x i ,Xi-2) 

= c{xi, Xi-2)b{xi, Xj_i)a(xj_ 2 , Xi-i) (19) 

Dividing the right hand side of equation 1191 by the expression in equation 
1171 we obtain 

9i-2,ifi,i-l + 9i-2,i-l9i— l,i = 9i-2,ifi-2,i-l 

So we end up with 

F\ = fi,i-lfi,i-2F-_ 2 + fi-l,i-2gi-l,%F?_2 + fi-2,i-igi-2,iF-_2 (20) 
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4.5 Rolling many times 



It is straightforward to iterate the above equation to bring the inverted arrow 
to the top. If the initial position is r, with rapidity x r , we obtain 



r-l 

E 

a=l 



r-l 

9a,r J^J fa.i 

i=l 



r — 1 \ r 

^=e 

k i=l 



a=l 



,/a,r i=1 



where we have used 



yr,r 
fr,r 



1 



Initially, we have the rows all labelled correctly, but the arrow on row r is 
inverted. 

After rolling the inverted arrow (in the remaining lattice) to the top row, 
as we did above, the top row (which has n — 1 vertices) is now 'frozen', in the 
sense that we know that all of its vertices are a vertices. Peeling these, we pick 
up a factor of 

N 
3=2 

where Xk is the rapidity that ended at the top. The partition function of the re- 
maining configurations can now be computed using Izcrgin's expression. Putting 
all contributions together, we obtain 



a=l 



N 

Yl a{x a -,Vj) 

U=2 



Ja.r 



i=l 



■>N-1 



where Zn-i [x <3l ,y 1 j , with underlined arguments, is Izcrgin's partition func- 
tion for an (N— 1) x (TV— 1) lattice, with the same assignment of horizontal and 
vertical rapidities, as those of Zn {y} j , but with x a and y\ missing. 

The result of Bogoliubov, Pronko and Zvonarev for the boundary 1- 
point function, H^, in their definition, can be produced by multiplying the left 
and the right 1-point functions computed above, and normalising the product 
by the partition function of the NxN model, to obtain 



H r N [{x},{y}} 



E 

a=l 



{{x},{y}) 



]J b ( x r,yi) 



N 



c(x r ,yi) 



' N 




r 




Y[a(x a ,yj) 




9a,r TT , 

f 11 a - i 

Ja.r i=1 


Zn-i 


.3=2 







[xa,^] 



(21) 
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4.6 Rolling down 



It is possible to repeat the above exercise by rolling downwards, rather than 
upwards. The resulting expression will look different from that obtained above. 
In particular, the peeled rows will be products of b rather than a weights. We 
leave it as an (easy) exercise to the reader to show that the resulting expression 
is trivially equivalent to the first. 

5 Simple boundary 2-point functions 

We are now ready to extend the above arguments to simple boundary 2-point 
functions. There are four cases to consider: 

1. The second inverted arrow is on the same boundary as the first. 

2. The second inverted arrow is on one of the two boundaries adjacent (rather 
than opposite) to that of the first. 

3. The second inverted arrow is on the opposite boundary to that of the first, 
but not directly opposite to it. 

4. The second inverted arrow is on the opposite boundary to that of the first, 
and directly opposite to it. 

Case 1 arises when an N xN domain wall lattice is sliced into two parts with 
one part having 2xN vertical (or horizontal lines) while the other part has the 
rest. 

Cases 2, 3 and 4 arise in more complicated dissections of the original domain 
wall lattice (for example, into more than 2 parts). We deal with them for 
completeness. 

5.1 Case 1 

A non-trivial example of a simple 2-point function, with both inverted arrows 
on the same side, can be drawn schematically as follows 




r> n 

— 

r> n 

— 

f V 



Figure 19 A simple boundary 2-point 
function. 
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If the initial positions of the inverted arrows are n and r 2 , where n < r 2 , 
with rapidities x ri and x r2 , respectively, we can roll the upper arrow to the top, 
then roll the lower arrow to the row below the first, and obtain 



ri,r 2 



E 



ai — 1 



/ai,n 



j J /ai,»i 



JJo(a; ai ,2/j) 

U=3 



pr 2 
r 2 



and 



^ - E 



9a 2 , 



fa 2 ,r 



n /- 



' AT 

J^a(x Q2 , yi ) 

J=3 



The full result is obtained by composing the above expressions. The extension 
of the above to (N-n)xN lattices with n inverted arrows, all on the same 
boundary, is clear. 



5.2 Case 2 

In this case, there are two types of configurations to consider. In both types, 
we are dealing initially with an NxN lattice. The first type can be drawn 
schematically as 



UN Vc V2 




N c 2 



Figure 20 A simple boundary 2-point function with inverted ar- 
rows on orthogonal boundaries. 

In the second type, we consider a horizontal inverted arrow pointing to the 
right boundary, and a vertical inverted arrow pointing away from the lower 
boundary. Both types can be treated analogously (leading to analogous results 
with minor differences). In the following we will deal with the first and leave 
the second as an exercise to the reader. 

If the positions of the inverted arrows are row r and column c, then we can 
write the result in two steps. First, we roll the inverted arrow on the right 
boundary all the way to the top to obtain 7 

7 In this case there is no simple frozen upper row that can be peeled. 
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This leaves us with the factor corresponding to the lattice 




N. 



N c 2 

Figure 21 The right inverted arrow has been rolled all the way to 
the top. 

We need to roll the vertical inverted arrow all the way to the left. This can 
be done using exactly the same method that was used above to roll a horizontal 
inverted arrow to the top. In this case, we need to 'inject' an a vertex into the 
lattice from below and use the Yang Baxter equation to thread it all the way 
to the top. The result is exactly the same as in the case of horizontal arrows 
given that the extra vertex will move through the lattice in the same direction 
as the orientations of both vertical rapidities that end up swapping positions. 
The resulting configurations 
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1 



V2 



N. 
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Figure 22 The upper inverted arrow has been rolled all the way to 
the left. 
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can be evaluated by inspection to be 



z?ca 



N 

£ 



N 



N-l 



\ \ a(x a ,yj) 



3 = 2 
3*0 



N 



U=2 



where we have defined 



- _ a{y ll y l+ i) 
o{yi,yi+i) 



5.3 Case 3 



(22) 



Zn-2 [x.i,y v Xa,y ) ( 23 ) 
c(yi,yi+i) 



This is a direct extension of case 1. We take the inverted arrows to be those on 
the right and left boundaries. 



r 2 o- 



N 



N 4 3 

Figure 23 A simple boundary 2-point function with inverted ar- 
rows on opposite boundaries, but not directly opposite to one an- 
other. 

Rolling the higher inverted arrow, then the lower one, we end up with the 
configuration 
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Figure 24 TTie inverted arrow have been rolled all the way to the 
top two rows. 



which can be evaluated by inspection to be 



where 



M = 1 



2=2 



ai—1 



9auri TT e 
f 11 Jai, 



' TV 

Y[a(x ai , yj ) 

J=3 



[rlF) 



U=3 



where we have used a self-explanatory notation to indicate the position and 
rapidity labels of the inverted arrow on the left boundary. 

5.4 Case 4 

In this case, we need to deal with inverted arrows that are directly opposite to 
each other. 
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Figure 25 A 2-point function with two inverted arrows on oppo- 
site boundaries, and directly opposite to one another. 

It is clear that in this case, the above rolling argument needs to be slightly 
modified, as there in no way to multiply by an a weight from one side, and 
thread that to the other side to roll an inverted arrow. 

We wish to show that this case can be reduced to a sum over cases of type 
3, discussed above. Consider two configurations drawn schematically as follows 
(all external arrows that are not drawn are meant to be the same in both 
configurations, and all internal arrows are summed over) 




Figure 26 The figure on the left is a schematic presentation of 
the set of configurations with directly opposite inverted arrows that 
we wish to roll up. The figure on the left is identical to the first, 
but the inverted arrow on the right is one row higher. 

The configuration on the left is the one that we are interested in. The 
configuration on the right is auxiliary. Multiply the left configuration by a b 
vertex, the right configuration by a c vertex and consider the sum of the two 
results 
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Figure 27 Multiplying the graph on the left above by a type b ver- 
tex, and that on the right by a type c vertex, and adding the results. 

In the sum, all internal arrows, in the external loop on the right, are now 
summed over, and the the sum can be represented in terms of the single graph 
on the left of equation |2HI Applying the Yang Baxter equation to weave that 
external vertex through the lattice we obtain the graph on the right of figure |2~51 




Figure 28 All arrow orientations that are allowed in the external 
loop on the right of the left figure are summed over, so we can use 
the Yang Baxter equation to weave that external vertex through the 
lattice from the right to the left. 

The resulting term on the right hand side of equation 1281 can now be decom- 
posed as 




Figure 29 Showing the content of the arrow configurations of the 
external loop on the left. 

Dividing all terms by b(xi,Xi-\), moving a term from the left hand side to 
the right, and using the antisymmetry of the b weights, we obtain 

[in) = [i-iFli] +9i,i-x [r'pr 1 ] +«-!,< (r 1 ^) (24) 

Equation El says that the problem with two opposing inverted arrow at row 
i can be reduced to the same problem at row (i — 1), plus two boundary 2-point 
functions of the type discussed in Case 3 above. 
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Using equation 1241 repeatedly, boundary 2-point functions with opposite in- 
verted arrows can be computed in terms of a sum boundary 2-point functions 
of type 3. 

Remark 4 One can show by direct calculation that \Ff vanishes, and hence 
[lF$) =92.1 [iFi] + gi , 2 [if?) 

5.5 Boundary n-point functions 

Simple boundary n-point functions can be constructed almost mechani- 
cally by repeated application of the above arguments. All that is needed is to 
develop conventions to handle the notational complexity of the resulting expres- 
sions. 

Composite boundary n-point functions will require summing over those 
inverted boundary arrows that are not frozen, which signals a proliferation of 
terms. However, there is hope that such summations can be simplified. 

In Bogoliubov et al evaluate not only boundary 1-point functions, but 
also the boundary spontaneous polarizations which are basically sums over many 
1-point functions. Bogoliubov et al use the algebraic Bethe ansatz to obtain 
expressions that do not contain more terms than a typical 1-point function 
evaluated at the same point as the boundary spontaneous polarization. 

This leads us to believe that sums that arise in computing 'composite' bound- 
ary n-point functions can also be performed, leading to relatively simple expres- 
sions. 

5.6 Bulk correlation functions 

It is clear how the bulk functions can be obtained as sums over products of 
boundary functions. In the simplest case of a 1-point function of a frozen arrow 
located deep in the lattice, the result is a product of two composite boundary 
correlation functions. 

It is highly unlikely that such (typically large) sums over products can be 
used to deduce physical properties of correlation functions. To get manageable 
results, we need to develop efficient methods to evaluate the composite boundary 
functions, which brings us back to the comments made in the above paragraphs. 

6 Epilogue 

The point of this work was to show that the result of Bogoliubov et al can re- 
produced using simple arguments. We also wished to show the same arguments 
apply, basically with no modification to higher simple boundary functions. The 
expressions that we obtain become cumbersome for higher n- functions, but they 
can be written basically by inspection of corresponding graph. 

We hope that there are ways to simplify these expressions by performing 
suitable summations. 
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